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In a recent work, Wyss, et.al. [Phys. Rev. Lett., 98, 238303 (2007)] have noted a property of 
'soft solids' under oscillatory shear, the so-called strain-rate frequency superposition (SRFS). We 
extend this study to the case of soft solids under large-amplitude oscillatory shear (LAOS). We show 
results from LAOS studies in a monodisperse hydrogel suspension, an aqueous gel, and a biopolymer 
suspension, and show that constant strain-rate frequency sweep measurements with soft solids can 
be superimposed onto master curves for higher harmonic moduli, with the same shift factors as 
for the linear viscoelastic moduli. We show that the behavior of higher harmonic moduli at low 
frequencies in constant strain-rate frequency sweep measurements is similar to that at large strain 
amplitudes in strain-amplitude sweep tests. We show surface plots of the harmonic moduli and the 
energy dissipation rate per unit volume in LAOS for soft solids, and show experimentally that the 
energy dissipated per unit volume depends on the first harmonic loss modulus alone, in both the 
linear and the nonlinear viscoelastic regime. 

PACS numbers: 83.60.Df, 83.80.Hj, 83.80.Kn 



I. INTRODUCTION 



Linear viscoelasticity[l[ is limited to small strains or 
strain rates, and assumes a constant dynamic viscosity 
and zero normal stress differences in shear flows. In 
order to probe their linear viscoelastic response, mate- 
rials are frequently subject to sma/Z-amplitude oscilla- 
tory shear, wherein an applied sinusoidal strain 7(i) — 
^Qsin{u}t) (here 70 is the strain amplitude, w is the an- 
gular frequency, t is the time) results in a stress response 
ait^Lo) = ^Q[G'i{uj)sin{ujt) + G'({io)cos{ujt)] at the same 
angular frequency as the input. Here G' and G" are the 
so-called 'storage' and 'loss' modulus L], and are func- 
tions of the angular frequency. By contrast, the resultant 
stress from a Zar(7e-amplitude oscillatory shear (LA0S)12l- 
[T]| test contains higher harmonics which may be inter- 
preted in terms of harmonic moduli G^ and G^ (subscript 
n refers to the n-th harmonic, see Section in(A) for the 
definition) , which are functions of both the strain ampli- 
tude and the angular frequency. The moduli G" may be 
physically interpreted in terms of the energy dissipated 
per unit volume per cycle of strain oscillation (see Sec- 
tion III(E) below). In a recent work, Ewoldt, etaL(Tl]| 
have assigned a physical meaning to the third harmonic 
moduli in terms of a deviation from linear viscoelastic 
behavior by decomposing the stress response into a su- 
perposition of an elastic and a viscous component which 
are then expanded in terms of Chebyshev polynomials of 
the first kind. The signs of the third-order Chebyshev 
coefficients have been interpreted in terms of strain stiff- 
ening/softening of the material using the elastic stress 
component, and in terms of shear thickening/thinning 
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using the viscous stress component. Meissner[125 car- 
ried out a comprehensive rheological study of three low- 
density polyethylene samples which were practically in- 
distinguishable through usual characterization methods, 
and in their linear viscoelastic response, but behaved dif- 
ferently when subject to non-linear deformation processes 
such as film blowing and extrudate swell. Apart from 
their obvious utility in industrial processes, LAOS stud- 
ies are of fundamental interest in rheology, and may lead 
to the development of more representative constitutive 
models. 

In a recent work, Wyss, etaZ.[l3], have shown a supris- 
ing feature in 'soft solids' (compressed emulsions, con- 
centrated suspensions, foams, pastes, gels, etc.), wherein 
constant strain-rate frequency sweep (SRFS) measure- 
ments of G'l and G" may be superimposed onto master 
curves, with the structural relaxation time t(7o) of the 
material showing an inverse power-law dependence on 
the strain-rate amplitude 70 (at high strain-rate ampli- 
tudes), viz. t(7'o) cx 7o~'' (j^ > is the exponent in the 
power-law). In typical soft solids, the linear viscoelastic 
moduli cross over at a frequency which is too low to be 
probed via conventional rheological tests, and the utility 
of the SRFS procedure is that by increasing the strain- 
rate amplitude, the crossover frequency can be shifted to 
a range where it may be directly probed. Apart from this, 
the SRFS procedure supplies a unified picture of the slow 
relaxation dynamics [l5| in a wide range of soft materials 
and is of intrinsic interest. Similar results have been con- 
firmed in a micelle- forming triblock copolyme r 116] , aque- 
ous foams [13], and dough mixed with gluten [ij. How- 
ever, the results in Wyss, et.aZ.fl^ have only been inter- 
preted in terms of linear viscoelastic moduli, although 
the authors have carried out several tests with strain 
amplitudes in the non-linear viscoelastic regime of their 
samples. 

In this paper, we present results from a systematic 
LAOS study of three representative soft solids, and show 
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that higher harmonic moduU in SRFS measurements can 
also be rescaled onto master curves, with the same shift 
factors as for the hnear viscoelastic moduh. We show 
that the higher harmonic moduh in SRFS measurements 
at low frequencies are similar to the moduli calculated at 
large strain amplitudes in strain-amplitude sweep tests. 
We show surfaces of the harmonic moduli and the energy 
dissipation rate per unit volume in LAOS and explicitly 
verify an earlier theoretical result Q that the energy dissi- 
pated per unit volume in a viscoelastic material depends 
on the first harmonic loss modulus alone in both the lin- 
ear and the nonlinear viscoelastic regime. In addition, 
we supply an unambiguous and detailed prescription for 
carrying out LAOS studies with soft solids. 



II. EXPERIMENT 

A. Instrumentation 

Our experiments were carried out at room temper- 
ature on the strain-controlled Advanced Rheometric 
Expansion System-2000 (ARES-2000, TA Instruments, 
United States') rheometer. We used a cone-plate mea- 
suring system[19j with a cone diameter of 25 mm (cone 
angle=0.1 rad). The rheometer permits acquisition of 
DC voltage signals from the torque transducer (for mea- 
surement of torque) and the optical encoder (for mea- 
surement of the motor angular deflection) through BNC 
connectors in the rear panel of the instrument. These 
unprocessed voltage signals (in the range ±5 volts) are 
not noise- filtered, or corrected for inertia and compli- 
ance of the torque transducer. Data was acquired at 16- 
bit analog input resolution through an analog-to-digital 
card (NT PCI-6014, National Instruments, United States) 
coupled with a Labview (National Instruments, United 
States) code at a sampling rate[2()i] of 10'^ points per sec- 
ond. A small DC offset was subtracted from the acquired 
oscillatory signal, and the signal was filtered for noise 
using a Savitzky-Golay filter[2l|. The signal was then 
calibrated [22] to find formulae which were used to con- 
vert the voltage values to quantities of physical interest. 
The calibration curves used were y = 0.046a; {x in volts, y 
in Nm) for the torque and y = Q.lx {x in volts, y in radi- 
ans) for the deflection angle. The values of the stress (in 
Pascals) and the strain were calculated from the torque 
and the deflection angle respectively using conversion fac- 
tors appropriate to the measuring-system geometry and 
torque transducer employed. 



B. Material synthesis and characterization 

The materials used in our tests were Poly N- 
isopropylacrylamide (PNIPAM), Xanthan gum, and 
Brylcreem Wetlook Gel. After loading the sample on 
to the rheometer, a thin layer of silicone oil (SF 1000, 



GE Bayer Silicones, India) was applied at the edges of 
the sample to prevent drying. 

PNIPAM synthesis was carried out in a double-jacket 
glass kettle reactor by a free-radical polymerization 
reaction j23||. The temperature was controlled by a water- 
circulating bath (Cool Tech 320, Thermo-Electron Corp., 
United States). 600 ml of distilled, deionised 'Milli- 
pore water' (Milli-Q Gradient AlO, Millipore, United 
States) purged with nitrogen gas for 1 hour, was used for 
preparing the solutions. The monomer NIPAM (Acros 
Chemicals, Belgium), the cross-linker N,N'-methylene- 
bisacrylamide BIS (Sigma- Aldrich, United States) and 
the initiator potassium peroxodisulfate KPS (Sigma- 
Aldrich, United States) were recrystallized from appro- 
priate solvents, and vacuum dried at room temperature 
for 4 hours. 7.87g of NIPAM, 0.39g of BIS, and 0.15g 
of the stabiliser sodium dodecyl sulfate (Merck, United 
States) were mixed in 48Gml of water at 25°C, and stirred 
using an overhead stirrer (RZR 2051control, Heidolph In- 
struments, Germany) at 100 rpm for 30 minutes under 
an inert atmosphere. 0.3g of KPS in 20 ml Millipore 
water, was added to the reaction mixture at 70°C, and 
stirred at 300 rpm, the reaction being allowed to pro- 
ceed for 4 hours. The temperature was then reduced to 
25°C, and the reaction mixture stirred overnight at 100 
rpm. Finally, the dispersion was dialysed (using dialysis 
bags having a molecular weight cutoff of 10000 g/mol) 
against Millipore water for two weeks. The dialysed sam- 
ple was frozen using liquid nitrogen, freeze-dried in a 
freeze drier (Heto Power Dry LL3000, Thermo-Electron 
Corp., United States) for 8 hours and stored in a dessi- 
cator. The prepared aqueous PNIPAM dispersion had 
a concentration of 14.09% w/w. The hydrodynamic ra- 
dius of PNIPAM was found to be 63.9 nm by Dynamic 
Light Scattering (DLS) on a Particle Size Analyser (BIG 
90Plus, Brookhaven Instrument Corp., United States). 

Xanthan gum (4% w/w) was prepared by dissolving 
4g of Xanthan gum (Sigma- Aldrich, United States) in 
lOOg Millipore water, under constant stirring using a Hei- 
dolph overhead stirrer for 8 hours. 7 mg of sodium azide 
(Merck, United States) was added to prevent bacterial 
growth. 

Brylcreem Wetlook Gel (Godrej Sara Lee Ltd., India) 
was purchased off-the-shelf and used as received. 



III. RESULTS AND DISCUSSION 

A. Large-amplitude oscillatory shear (LAOS) 

Using results from continuum mechanics, 
Christensen[23| has shown that for an imposed si- 
nusoidal strain ^[t) = josin tot (here 70 is of arbitrary 
magnitude), the resultant stress in an isotropic material 
may be written as an odd-harmonic Fourier series: 

cr(t; w, 70) = 70 E„=i,3,5,... [(^"(w, 7o)sm nujt 

+G'l^{uj,-1o)cos nut], (1) 
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where and G" are n-th harmonic moduh which are, 
in general, functions of both the angular frequency and 
the strain amplitude. 

We define the n-th harmonic moduli as 

G'n{^,lo) = —COS $„, 

71 

G::(c.,7o) = ^szn$„, (2) 
71 

where $„(cli,7o) = 4>n.a ~ n(j)^ are the phase angles for 
the n-th harmonic moduli o- and 0^ are the phase 
angles of the n-th harmonic in the stress and the first 
harmonic in the strain Fourier series, respectively). Here 
cr„(u;,7o) is the n-th harmonic in the stress amplitude 
spectrum, and 71 is the first harmonic in the strain am- 
plitude spectrum 25]. Note that in our definition of the 
moduli [Eqs. ([H] we distinguish between 70, the re- 
quested strain amplitude, and 71, the measured ampli- 
tude of the first harmonic in the strain amplitude spec- 
trum - these quantities differ by a small amount, which 
was empirically found to vary with the torque transducer 
compliance. 



B. Data analysis 

The general stress response in LAOS [Eq. ([1])] may be 
rewritten as cr(i; oj, 70) = X)n=i 3 5 d'nsin{nujt + 4>n,a), 
where cr„ and 4>n,<T are both real-valued functions of lo 
and 7o. In LAOS theory, the strain is assumed to be a 
sine wave with zero phase, but an experimental strain sig- 
nal ['y(t) = 7osm(wi-|-(/).y)] has a non-zero phase angle (pj 
whose value depends on the instant of time when data ac- 
quisition commences. Hence, it is useful to apply the fol- 
lowing transformation: t = t' — (p-y/uj, and obtain ^{t') — 
josin{ujt') and a{t';uj,^o) — 3 5 ansin{nujt' + 

4'n,a — n(j)-y)- Therefore, data analysis must be carried 
out on complete oscillation cycles of the strain signal, 
with zero phase. 

Data analysis was carried out using Matlab (The 
MathWorks Inc., United States). The full- width at half 
maximum (FWHM) of the sine function[26|] is approxi- 
mately 1/(2T)[2^, where T is the duration of the time 
series. We choose the duration of our runs to correspond 
to a value of FWHM equal to 0.5% of the imposed angu- 
lar frequency, this corresponds to a run with at least 50/7r 
cycles of oscillations. In the ARES-2000 rheometer, the 
measured strain amplitude in oscillatory shear, reaches 
the requested value 70 after a few transient cycles of os- 
cillation. We ensure that our signal is processed from 
the instant of time when the measured strain amplitude 
is close to the requested value, and use a zero-crossing al- 
gorithm to extract 50 cycles of oscillation from the strain 
and the stress signal (which have the same time-base). 
Note that for the materials under study, the amplitude 
of the measured oscillatory torque signal does not decay 
appreciably during the course of the experiments. We 



apply a discrete Fourier transform on the extracted sig- 
nals, and use only the first half of the Fourier-transformed 
dataset (remainder of the dataset is redundant for real 
input data). Although in theory, the higher harmonic 
frequencies are integer multiples of the fundamental fre- 
quency, in an experimental stress signal the measured 
harmonic frequencies located through peak values in the 
stress amplitude spectrum, differ by a small amount from 
the integer multiples. Hence, we use a code which locates 
peak values in the stress amplitude spectrum and explic- 
itly finds the frequencies corresponding to the harmonics. 
Finally, we calculate the phase angles and thereafter the 
harmonic moduli of interest. 

In this study, we do not apply oversampling [28j tech- 
niques to improve the signal-to-noise ratio, and the res- 
olution of the moduli beyond the fifth was found to be 
poor, due to noise and the low sampling rate. In ad- 
dition, although the first harmonic moduli reported by 
the proprietary software accompanying the ARES-2000 
rheometer are reliable up to a; = 100 rad/s (the up- 
per limit for the machine), we found that the calculated 
higher harmonic moduli are of relatively low resolution 
beyond w « 50 rad/s. Therefore, we restrict our study of 
the higher harmonic moduli up to the fifth order, with a 
maximum angular frequency of 50 rad/s. It is important 
to note that there is no automated means of obtaining 
higher harmonic moduli from the ARES-2000 rheometer, 
therefore, each point in our graphs (with the exception of 
Figs. [T]and[2]) is calculated from an independent oscilla- 
tory shear test with a fixed angular frequency and strain 
amplitude. 



C. Strain-rate frequency superposition (SRFS) in 
LAOS 

Most of the results that follow are from experiments 
carried out with our highly monodisperse PNIPAM sam- 
ple, although a few results are also reported from studies 
with Xanthan gum. The results from tests carried out 
with Brylcreem are reported in one instance. 

In Fig. [IJa), we plot the first harmonic moduli G'l 
and G" as a function of varying strain amplitude 70 with 
the angular frequency w = 1 rad/s, using PNIPAM. The 
plot shows features characteristic of soft solids, with a 
pronounced peak in G'{ at intermediate frequencies and 
a corresponding fall in G'l , indicative of a breakdown of 
internal structure, and resultant liquid-like behavior at 
higher strain amplitudes Ts*] . G'{ is a measure of the en- 
ergy dissipated per unit volume per cycle of strain oscil- 
lation (see Section III(E) below) and the growth towards 
a peak value may be interpreted in terms of energy loss 
(as heat), on account of the work done to break inter- 
nal structures in the material. In Fig. [IJb), we plot the 
first harmonic moduli as a function of uj for 70 = 0.01, a 
strain amplitude within the linear viscoelastic regime of 
PNIPAM. In the range of frequencies probed, the sample 
shows solid-like behavior with G[ greater than G'/, but 
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FIG. 1: (a) Plot of the first harmonic moduli G'l (solid circles) 
and Gi (open circles) as a function of varying strain ampli- 
tude 70 with Lj = 1 rad/s using PNIPAM. 
(b) Plot of the first harmonic moduli G'l (solid circles) and 
G'l (open circles) as a function of varying angular frequency 
oj with 70 — 0.01 using PNIPAM. The dashed line is propor- 
tional to Li;*^"^* (see the text for discussion). 



the small upturn in G" at the smallest frequencies hints 
at the possibility of a crossover of the moduli at a much 
lower frequency. The dashed-line curve in this figure is 
discussed below. 

In Fig. EJa), we plot a few cycles of sinusoidal oscilla- 
tions of the strain signal as a function of the time from a 
LAOS test at 70 = 1.2, w = 1 rad/s using PNIPAM. In 
Fig. El^b), we plot the resultant stress signal as a function 
of the time. Note the non-sinusoidal nature of the stress 
signal, indicative of higher harmonics in the stress am- 
plitude spectrum. In Fig. ^c), we plot the strain power 
spectrum as a function of ui for 50 oscillation cycles 
of the strain signal, a few oscillation cycles of which are 
shown in Fig. [DJa). Here Ph(w) = 2|i?(w)p is the one- 
sided power spectrum of the real- valued time-series h(t) 
with Fourier coefficients H{uj). The higher harmonics in 
the measured strain power spectrum may be considered 
as part of broadband noise. Our definition of the moduli 
[Eqs. ([2])] and resulting calculations do not account for 
higher harmonics in the strain signal. In Fig. [H^d), we 
plot the response stress power spectrum as a func- 
tion of oj. Prominent odd harmonics [2^ can be seen at 
frequencies which are close to integer multiples of the ap- 
plied oscillation frequency uj = 1 rad/s along with much 
smaller even harmonics, a result first noted by Krieger 
and LiuQ. Even harmonics in the stress power spec- 



trum have been hypothesized [30j to occur in oscillatory 
shear flows showing wall-slip, or due to secondary flows 
in the plate gap. In our experiments, the ratio of the sec- 
ond to the flrst harmonic in the stress power spectrum 
is of order 10~® or smaller, even harmonics are therefore 
neglected. 

In Fig. [Sl^a), we plot the first harmonic moduli from 
SRFS tests with values of the strain-rate amplitude 'Jq ~ 
7oa; in the range [0.9 s~^,4.2 s~^], as a function of uj us- 
ing PNIPAM. Note that the moduli crossover frequency 
is shifted in the direction of increasing w, with larger val- 
ues of T'o and the shapes of the moduli curves are similar. 
In Fig. [3jb), we plot the rescaled moduli as a function 
of the rescaled angular frequency. The values of the shift 
factors ai(7'o) and 61 (t'o) are listed in Table I and the 
moduli curves for the highest value of T'o are taken as the 
reference curves for the purpose of shifting. Note that al- 
though G"/ is strictly positive (see Section III(E) for dis- 
cussion), there is no restriction on the sign of the other 
moduli, which is why we plot all our moduli curves on a 
semilogarithmic scale for consistency. The plot confirms 
the validity of the SRFS procedure for linear viscoelastic 
moduli, as reported in Wyss, et.aL[l^. To ensure that 
the recorded torque signal was well within the capabili- 
ties of the machine, our experiments were not carried out 
at small strain-rate amplitudes, which is why we were un- 
able to confirm the low-strain rate dependence of 6(70)1 
asymptotically expected to approach unity [11]. A discus- 
sion regarding the high-frequency response of the moduli 
may be found at the end of this Section. 

Our results are, however, not restricted to the linear 
viscoelastic regime, and in Figs. Ilja),(b), we show the 
corresponding rescaled third and fifth harmonic mod- 
uli as a function of the rescaled angular frequency from 
the same test as reported in Fig. [3] The values of 
the shift factors are listed in Table I. Note the strik- 
ing agreement between the values of the vertical shift 
factors (03(70)1 o-bilo)) with ai(7'o) and the horizontal 
shift factors (63(70), 65(70)) with 61(70). In Fig. IDJc), 
we plot the shift factors as a function of 7*0. The shift 
factors a„(7'o) (n = 1,3,5) were found to be of order 
unity, while the shift factors bnijo) [i^ — Ii3j5) showed 
a power-law dependence on 7-0 (6„(7o) oc 7*0 '^), with an 
exponent 1/ = 0.89 ± 0.01. The shift factor 61(70) is 
postulated[l^ to depend inversely on the structural re- 
laxation time, viz. 61(70) oc 1/t(7'o), and our results for 
bniio) = 1,3,5) are in agreement with the exponent 
v fa 0.9 reported in Wyss, et.aZ.[l^. Physically, this sug- 
gests that soft solids have internal structures that relax 
faster in the presence of an imposed strain-rate (the in- 
verse of which essentially "sets" the structural relaxation 
time of the material). 

Analogous master curves from tests using Xanthan 
gum are shown in Fig. [5] for the first, third and the 
fifth harmonic moduli. The shift factors for the higher 
harmonic moduli were again found to agree (see Table 
I) with the shift factors for the first harmonic mod- 
uli. The horizontal shift factors 6„(7'o) (n = 1,3,5) 
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FIG. 2: (a) Plot of a typical strain signal 7(f) as a function of 
the time t from an oscillatory shear test with 70 = 1.2, cj = 1 
rad/s using PNIPAM. 

(b) Plot of the resultant stress signal a{t) as a function of the 
time t with parameter values as in (a). 

(c) Plot of the power spectrum (see the text for the defini- 
tion) of the strain signal (a few oscillations are shown in (a)) 
as a function of the angular frequency ui. 

(d) Plot of the power spectrum Pc of the stress signal (a few 
oscillations are shown in (b)) as a function of the angular 
frequency u. 



FIG. 3: (color online) (a) Plot of the first harmonic moduli G'l 
(solid circles) and G" (open circles) as a function of the an- 
gular frequency uj from constant strain-rate frequency sweep 
measurements at strain-rate amplitudes 70 = 4.2 (blue), 
2.1 s"^ (red), 1.2 s"^ (green), 0.9 s"^ (pink) using PNIPAM. 
(b) Plot of the scaled first harmonic moduli G''i/ai(7o) (solid 
circles) and Gi/ai(7o) (open circles) as a function of the 
scaled angular frequency 01/61(70) shifted onto a single master 
curve. See Table I for the values of 01(70) and 61(70). 



showed a power-law dependence on 7-0 with an exponent 
1/ = 0.89 ±0.05. Results from our Brycreem tests are not 
shown in this work, apart from the shift factors listed in 
Table I, and serve to confirm our primary result. In the 
case of Brycreem, the power-law dependence of bn{io) on 
70 followed, with = 0.99 ± 0.11. 

In Fig. [Il^b), the first harmonic moduli cross over at 
a frequency which is not easily accessible, however, the 
shape of the SRFS master curve in Fig. ^h), suggests 
that the low frequency behavior of the moduli is simi- 
lar to that at high strain amplitudes (compare with Fig. 
[IJa)). We confirm an analogous result for higher har- 
monic moduli: Figure IHl^a) is a plot of the third har- 
monic moduli G3 and G3 as a function of uj from an 
SRFS test with = 1.2 s"^ using Xanthan gum, while 
Fig. is a plot from a varying strain amplitude test at 
a; = 5 rad/s. The plots are a remarkable demonstration 
of the 'reversed' nature of the higher harmonic moduli as 
reported from SRFS tests and from a strain-amplitude 
sweep test. Similar results have been confirmed for the 
fifth harmonic moduli. 

In Wyss, et.aL[l3|. it is shown that at high frequen- 
cies, the curves for G" for different strain-rate ampli- 
tudes can be rescaled onto a single master curve if a y/uj 
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FIG. 4: (color online) (a) Plot of the scaled third harmonic 
moduli 03/03(70) (solid circles) and 03/13(70) (open circles) 
as a function of the scaled angular frequency u)/bz{')o) shifted 
onto a single master curve from constant strain-rate frequency 
sweep measurements at strain- rate amplitudes 7*0 = 4.2 
(blue), 2.1 (red), 1.2 (green), 0.9 s"^ (pink) using 
PNIPAM. See Table I for the values of 03(70) and 63(70). 

(b) Corresponding plot of 05/05(70) and 05/05(70) as a func- 
tion of Lu/b5{'yo). See Table I for the values of 05(70) and 
65(70). 

(c) Plot of the shift factors an(7'o) (solid squares) and 6n(7'o) 
(solid circles) as a function of the strain-rate amplitude 7-0 for 
the moduli G'n and Gn [n = 1 (blue), 3 (red), 5 (green)]. 



FIG. 5: (color online) (a) Plot of the scaled first harmonic 
moduli O'i/oi(7o) (solid circles) and Oi/oi(7o) (open circles) 
as a function of the scaled angular frequency u)/bi{'yo) shifted 
onto a single master curve from constant strain-rate frequency 
sweep measurements at strain-rate amplitudes 7*0 = 4.2 
(blue), 2.1 (red), 1.2 (green), 0.3 (pink) using 
Xanthan gum. See Table I for the values of 01(70) and 61(70). 

(b) Corresponding plot of 03/03(70) and 03/03(70) as a func- 
tion of aj/63(7o). See Table I for the values of 03(70) and 
63(70). 

(c) Corresponding plot of 05/05(70) and 05/05(70) as a func- 
tion of uj/bs{'yo). See Table I for the values of 05(70) and 
65(70). 



component is subtracted, the high-frequency linear vis- 
coelastic response being attributed to viscous flow along 
randomly oriented slip planes 'si']. However, in the an- 
gular frequency sweep for PNIPAM (see Fig. lUb)), the 
dashed-Hne fit to the high-frequency portion was found 
to be proportional to 0;" °^^. Clearly, our PNIPAM sam- 
ple has not reached the asymptotic high-frequency regime 
at w = 100 rad/s (the upper-limit for ARES-2000) where 
the scaling law applies. For Xanthan gum and Brylcreem, 



similar results were obtained, with a high-frequency scal- 
ing for G'{ proportional to uj^-'-^^^ and a;°'^® respectively. 
Therefore, wc did not attempt to correct for the observed 
lack of collapse of G'{ onto a master curve in Figs. [3Jb) 
and [5l^a) at high frequencies. Apart from this, an ad- 
ditional effect which must be accounted for is the in- 
ertia of the torque transducer and the variation of the 
torque transducer compliance due to vibrations trans- 
mitted through the frame of the rheometer, under high- 
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TABLE I: Vertical {a„{"fo)) and horizontal (6,1(70)) shift fac- 
tors for the harmonic moduli G'„ and G'n {n = 1,3,5) for 
different strain-rate amplitudes 70 using PNIPAM, Xanthan 
gum and Brylcreem. 
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0.83 
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0.05 
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0.12 
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frequency oscillations of the motor |32l|. The ARES-2000 
rheometer incorporates a hardware correction for these 
effects and there is no means by which the corrections 
may be implemented through an external software. The 
shift factors listed in Table I, however, remain unaffected. 

D. Surfaces of harmonic moduli 

Most LAOS studies to-date have been restricted to the 
study of stress-strain curves 0, H, H, [ll| or the ratios of 
the harmonics of the stress amplitude spectrum [5, 0, 
In occasional studies the harmonic moduli have been 
calculated in a strain amplitude or an angular frequency 
sweep test. However, the complete form of the moduli is 
seen clearly in a surface plot, obtained from oscillation 
tests at different values of 70 and w. In an early study, 
Reimers and Dealy|j9|] used a sliding-plate rheometer and 
plotted surfaces of the third and the fifth harmonic mod- 
uli from LAOS studies with polystyrene. The plots in 
their study are qualitative at best and of low resolution. 

In Figs. [3 [H and[9l we show surface plots of the har- 
monic moduli. These plots were obtained from 256 inde- 
pendent oscillatory shear tests using Xanthan gum, with 
the values of 70 and uj spaced logarithmically on a 16 x 16 
grid, the elements along one diagonal having a strain- 
rate amplitude 7*0 = Jqlo — 2.1 s^^ (a value chosen from 
Table I). Although the first harmonic loss modulus G" 
is strictly positive (see Section III(E) below), the other 
moduli can be of any sign, which is why our surface plots 
are shown on a logarithmic scale in the (cj,7o) plane, 
and with a linear scale on the third axis. Note that the 
SRFS results in Wyss, ei.aZ.flJ] pertain to the study of 
curves of intersection of the hyperbolic sheets 70W = F (F 
are the chosen values of the strain-rate amplitudes, here 
F = 2.1 s"^) with the corresponding harmonic surfaces 
for G"]^(a;,7o) and G'({uj,^q). Similarly, the intersection 




FIG. 6: (a) Plot of the third harmonic moduli G'^ (solid cir- 
cles) and G3 (open circles) as a function of the angular fre- 
quency UJ with constant strain-rate amplitude 70 = 1.2 , 
using Xanthan gum. 

(b) Plot of the third harmonic moduli G3 (solid circles) and 
Gg (open circles) as a function of the strain amplitude 70 for 
angular frequency cj = 5 rad/s, using Xanthan gum. 

of the hyperbolic sheet F — 2.1 s^^ with the harmonic 
surfaces plotted in Figs. |8] and [9] would supply the curves 
marked with red circles in Figs. [5jb) and (c) respectively. 

We observe that the surface for G'{ has a local max- 
ima at (a;,7o) = (2.31 rad/s, 0.17), while the surfaces for 
G3, G3, G5, and Gg have local minimas at (w,7o) = 
(1.32 rad/s, 0.097), (1.14 rad/s, 0.17), (1.32 rad/s, 0.225), 
and (1.14 rad/s, 0.393) respectively. The moduli G" and 
G3 share the same value of 70, while the pairs (G3, G5) 
and (G3 , G5 ) share the same values of uj. This interesting 
coupling between the moduli requires further investiga- 
tion in surface plots with much higher resolution, and 
may be of interest to rheologists modelling soft glassy 
materials. 



E. Energy dissipation and dissipation rate 

The energy dissipated per unit volume per cycle of 
strain oscillation is e(i; 0^,79) = ^^'^^'^ ad'y^. Ganeri- 
wala and RotzQ have shown that on substituting the 
one-dimensional Green-Rivlin constitutive eauation[33j 
for the stress into this formula and assuming a sinusoidal 
strain one obtains e = ttjqG'Kuj, 70), an expression which 
is true for arbitrary values of 70. The second law of ther- 
modynamics requires that e > 0, i.e. G'{ is strictly pos- 



8 




FIG. 7: (color online) (a) Surface plot of the first harmonic 
modulus Gi(w,7o) as a function of strain amplitude 70 and 
angular frequency cj, using Xanthan gum. The color bar in- 
dicates magnitude of the modulus in units of Pascal, 
(b) Surface plot of the first harmonic modulus G'/ (01,70) as 
a function of strain amplitude 70 and angular frequency uj, 
using Xanthan gum. 



FIG. 8: (color online) (a) Surface plot of the third harmonic 
modulus G3(aj,7o) as a function of strain amplitude 70 and 
angular frequency uj, using Xanthan gum. The color bar in- 
dicates magnitude of the modulus in units of Pascal, 
(b) Surface plot of the third harmonic modulus 03(0;, 70) as 
a function of strain amplitude 70 and angular frequency ui, 
using Xanthan gum. 



itive, but places no restriction on the sign of the other 
moduli. It is plausible that the other harmonic moduli 
are involved in reversible exchanges of energy in LAOS. 

In Fig. [ini we plot the ratio e/{njQG'-l) (e being the 
area bounded by the stress-strain curve |34]), as a func- 
tion of the number of the point p (each point representing 
an independent test) for three different tests: A constant 
strain-rate frequency sweep test at = 0.9 {uj in 
the range [1,25] rad/s, 70 in the range [0.036,0.9]) us- 
ing PNIPAM, a frequency sweep test at 70 — 0.52 {uj 
in the range [1,25] rad/s) using Xanthan gum, and a 
strain-amplitude sweep test at w = 9.35 rad/s (70 in the 
range [0.042,2.1]) using Xanthan gum. In all three test 
configurations, we note that the points cluster about the 
expected value of 1, with a maximum deviation of ap- 
proximately 2%, the residual difference being attributed 
to experimental uncertainty in the measurement of G'(. 

The energy-dissipation rate per unit volume in LAOS 
can be shownfl^ to equal e(t;a;,7o) = aj7oG"(a;, 7o)/2. 
In Fig. [TIT a). we show the surface plot of the energy- 
dissipation rate per unit volume for Xanthan gum. It is 
interesting to note that the e surface grows monotonically 
with 7o and despite the material showing a solid-like 
response at high frequencies. On account of the logarith- 
mic scaling used for the axes and the orientation of the 
surface in Fig. Illf a). it may not be directly apparent 



that e is in fact a decreasing function at large uj along the 
curve of intersection of the surface t'o — uj^o — 2.1 
with the surface for e, as is shown in Fig. [TTJb). We may 
rewrite the expression for the energy dissipation rate per 
unit volume as e = 7p^G"/(2a;). In linear viscoelasticity, 
G'( (X UJ at small a;[l|, implying that e is constant along 
the SRFS curve (with constant T'o) at small uj, which 
accords with Fig. [TIT b) for T'o = 2.1 s~^. The dashed- 
line fit to the high-frequency (corresponding to low strain 
amplitude) portion of the e curve was found to be propor- 
tional to uj~'^-^^. For comparison, in linear viscoelasticity, 
the Maxwell model (see Ref. [l], p. 57) predicts e c>c l/w^ 
along the SRFS curve at high frequencies. 



IV. SUMMARY 

To conclude, we have presented results from a system- 
atic experimental study of soft solids under LAOS with 
special attention to SRFS measurements. Our results 
show the general applicability of the SRFS result for soft 
solids, specifically, that the SRFS curves for higher har- 
monic moduli can be superimposed onto master curves, 
with the same shift factors as for the linear viscoelastic 
moduli. We have also shown surface plots of the mod- 
uli and the energy-dissipation rate per unit volume in 
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FIG. 9: (color online) (a) Surface plot of the fifth harmonic 
modulus G5(a;,7o) as a function of strain amplitude 70 and 
angular frequency cj, using Xanthan gum. The color bar in- 
dicates magnitude of the modulus in units of Pascal, 
(b) Surface plot of the fifth harmonic modulus G'^{ijj,^o) as 
a function of strain amplitude 70 and angular frequency to, 
using Xanthan gum. 
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FIG. 10: (color online) Plot of t/ {ifyQC'l) as a function of 
the number of the point p (each point representing an inde- 
pendent test) from three different tests, a constant strain-rate 
frequency sweep test (red) at 7-0 = 0.9 using PNIPAM, 
a strain-amplitude sweep test (blue) a.t uj = 9.35 rad/s and a 
frequency sweep test (green) at 70 = 0.52, both using Xan- 
than gum. 



LAOS. We have shown that the energy dissipated per 
unit volume in oscillatory shear is governed by the first 
harmonic loss modulus alone. We hope that the results 
in this paper motivate further studies of LAOS in poly- 
mer melts, electro- and magnetorheological fluids, among 
other complex materials. 




FIG. 11: (a) (color online) Surface plot of the energy- 
dissipation rate per unit volume e(t; 01,70) = ctJ7oGi(aj, 7o)/2 
as a function of the angular frequency ui and the strain am- 
plitude 70, using Xanthan gum. The color bar indicates the 
magnitude of e in units of Pa/s. 

(b) Curve of intersection of the surface 70 = W)o = 2.1 
with the surface plot in (a). The dashed line is proportional 

to 
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